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Abstract 

Modelling multiple partonic interactions in hadronic events is vital for un- 
derstanding minimum-bias physics, as well as the underlying event of hard 
processes. A brief overview of the current Pythia 8 multiple interactions 
(MI) model is given, before looking at two additional effects which can be in- 
cluded in the MI framework. With rescattering, a previously scattered parton 
is allowed to take part in another subsequent scattering, while with enhanced 
screening, the effects of varying initial-state fluctuations are modelled. 
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1 Introduction 



The run-up to the start of the LHC has led to a greatly increased interest in the physics of 
multiple parton interactions in hadronic collisions. Existing models are used to try to get 
an insight into what can be expected at new experiments, extrapolating fits to Tevatron 
and other data to LHC energies pQ. Such extrapolations, however, come with a high level 
of uncertainty; within many models are parameters which scale with an uncertain energy 
dependence. There is, therefore, also the exciting prospect of new data, with which to 
further constrain and improve models. 

In terms of theoretical understanding, MI is one of the least well understood areas. 
While current models, after tuning, are able to describe many distributions very well, there 
are still many others which are not fully described. This is a clear sign that new physical 
effects need to be modelled and it is therefore not enough to "sit still" while waiting for 
new data. It is with this in mind that we look at two new ideas in the context of MI and 
their potential effects. 

With rescattering, an already scattered parton is able to undergo another subsequent 
scattering. Although, in general, such rescatterings may be relatively soft, even when 
compared to normal 2 — > 2 MI scatterings, they can lead to non-trivial colour flows which 
change the structure of events. Another idea is to consider partonic fluctuations in the 
incoming hadrons before collision. In such a picture, it is possible to get varying amounts 
of colour screening on an event-by-event basis. The question then is, what effects such new 
ideas would have on multiple interactions and how can they be included in the Pythia 
framework? 

In Section El a brief introduction to the existing MI model in Pythia 8 is given. For 
more comprehensive details about what is contained in the model, readers are directed 
to [2] and the references therein. In Sections [3] and HI an initial look at rescattering and 
enhanced screening is given. A summary and outlook is given in Section [5l 

2 Multiple Interactions in Pythia 8 

The MI model in Pythia 8 [3J is a model for non-diffractive events. It is an evolution of the 
model introduced in Pythia 6.3 |2j, which in turn is based on the model developed in earlier 
versions of Pythia. The earliest model jl] was built around the virtuality-ordered parton 
showers available at the time and introduced many key features which are still present in 
the later models, such as p± ordering, perturbative QCD cross sections dampened at small 
p±, a variable impact parameter, PDF rescaling, and colour reconnection. 

The next-generation model E] was developed after the introduction of transverse- 
momentum-ordered showers, opening the way to have a common p± evolution scale for 
initial-state radiation (ISR), final-state radiation (FSR) and MI emissions. The second key 
ingredient was the addition of junction fragmentation to the Lund String hadronisation 
model, allowing the handling of arbitrarily complicated beam remnants. This permitted 
the MI framework to be updated to include a more complete set of QCD 2^2 processes, 
with the inclusion of flavour effects in the PDF rescaling. 

The Pythia 8 MI framework also contains additional new features which are not found 
in previous versions, such as 
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• a richer mix of underlying-event processes (7, J/tp, Drell-Yan, etc.), 

• the possibility to select two hard interactions in the same event, and 

• the possibility to use one PDF set for hard processes and another for other subsequent 
interactions. 



2.1 Interleaved p± Ordering 

Starting in Pythia 6.3, ISR and MI were interleaved with a common p± evolution scale. 
In Pythia 8, this is taken a step further, with FSR now also fully interleaved. The overall 
probability for the i th interaction or shower branching to take place at p± = p^ is given by 

dV ( dV M \ , d^iSR , d ^FSR S 



dp± I d P J 



+ E 




, rp±i-i dV M1 dV lSR dV FSR \ , \ 

with contributions from MI, ISR and FSR unitarised by a Sudakov-like exponential factor. 
If we now focus on just the MI contribution, the probability for an interaction is given 

dp±i a nd dp± \ J P± o- n ddp' ± J 

where da/dp± is given by the perturbative QCD 2^2 cross section. This cross section 
is dominated by t-channel gluon exchange, and diverges roughly as dp\/p\. To avoid this 
divergence, the idea of colour screening is introduced. The concept of a perturbative cross 
section is based on the assumption of free incoming states, which is not the case when 
partons are confined in colour-singlet hadrons. One therefore expects a colour charge to be 
screened by the presence of nearby anti-charges; that is, if the typical charge separation is 
d, gluons with a transverse wavelength ~ l/p± > d are no longer able to resolve charges 
individually, leading to a reduced effective coupling. This is introduced by reweighting the 
interaction cross section such that it is regularised according to 

d ^ a|(pj) a|(pio+Pi) /on 

dpi pi fc + pi) 2 ' [) 

where p± Q (related to 1/d above) is now a free parameter in the model. 



2.2 Impact Parameter 

Up to this point, all parton-parton interactions have been assumed to be independent, such 
that the probability to have n interactions in an event, V n , is given by Poissonian statistics. 
This picture is now changed, first by requiring that there is at least one interaction, such 
that we have a physical event, and second by including an impact parameter, b. The default 
matter distribution in Pythia is a double Gaussian 

p(r) "^T exp (-^) + | exp (-^|)' (4) 
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such that a fraction (3 of the matter is contained in a radius ct 2 , which in turn is embedded in 
a radius a± containing the rest of the matter. The time-integrated overlap of the incoming 
hadrons during collision is given by 



after a suitable scale transformation to compensate for the boosted nature of the incoming 
hadrons. 

Such an impact parameter picture has central collisions being generally more active, with 
an average activity at a given impact parameter being proportional to the overlap, 0(b). 
While requiring at least one interaction results in V n being narrower than Poissonian, when 
the impact parameter dependence is added, the overall effect is that V n is broader than 
Poissonian. The addition of an impact parameter also leads to a good description of the 
"Pedestal Effect" , where events with a hard scale have a tendency to have more underlying 
activity; this is as central collisions have a higher chance both of a hard interaction and of 
more underlying activity. This centrality effect naturally saturates at p± hard ~ 10 GeV. 



In the original model, PDFs were rescaled only such that overall momentum was conserved. 
This was done by evaluating PDFs at a modified x value 



where the subscript i refers to the current interaction and the sum runs over all previous 
interactions. The original model was affected by a technical limitation in fragmentation; it 
was only possible to take one valence quark from an incoming hadron. This meant that the 
MI framework was limited to qq and gg final states and that it was not possible to have 
ISR from secondary scatterings. By introducing junction fragmentation, where a central 
junction is connected to three quarks and carries baryon number, these limitations were 
removed. This allowed the next-generation model to include a more complete set of MI 
processes and flavour effects in PDF rescaling. 

ISR, FSR and MI can all lead to changes in the incoming PDFs. In the case of FSR, 
a colour dipole can stretch from a radiating parton to a beam remnant, leading to (a 
modest amount of) momentum shuffling between the beam and the parton. Both ISR 
and MI can result in large x values being taken from the beams, as well as leading to 
flavour changes in the PDFs. If a valence quark is taken from one of the incoming hadrons, 
the valence PDF is rescaled to the remaining number. If, instead, a sea quark (q s ) is 
taken from a hadron, an anti-sea companion quark (q c ) is left behind. The x distribution 
for this companion quark is generated from a perturbative ansatz, where the sea/anti-sea 
quarks are assumed to have come from a gluon splitting, g — > q s q c . Subsequent perturbative 
evolution of the q c distribution is neglected. Finally, there is the issue of overall momentum 
conservation. If a valence quark is removed from a PDF, momentum must be put back 
in, while if a companion quark is added, momentum must be taken from the PDF. This 
is done by allowing the normalisation of the sea and gluon PDFs to fluctuate such that 
overall momentum is conserved. 




(5) 



2.3 PDF Rescaling 
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2.4 Beam Remnants, Primordial k± and Colour Reconnection 



When the p± evolution has come to an end, the beam remnant will consist of the remaining 
valence content of the incoming hadrons as well as any companion quarks. These remnants 
must carry the remaining fraction of longitudinal momentum. PYTHIA will pick x values 
for each component of the beam remnants, according to distributions such that the valence 
content is "harder" and will carry away more momentum. In the rare case that there 
is no remaining quark content in a beam, a gluon is assigned to take all the remaining 
momentum. 

The event is then modified to add primordial k±. Partons are expected to have a non- 
zero k± value just from Fermi motion within the incoming hadrons. A rough estimate 
based on the size of the proton gives a value of ~ 0.3 GeV, but when comparing to data, 
for instance the p± distribution of Z° at CDF, a value of ~ 2 GeV appears to be needed. 
The current solution is to decide a k± value for each initiator parton taken from a hadron 
based on a Gaussian whose width is generated according to an interpolation 

a{Q) = max ^a min , 1 + g^g j > ( 7 ) 

where Q is the hardness of a sub-collision, cr m in is a minimal value (~ 0.3 GeV), is a 
maximal value that is approached asymptotically and Qi is the Q value at which a(Q) is 
equal to half a^. The recoil is shared among all initiator and remnant partons from the 
incoming hadrons, and the k± given to all daughter partons through a Lorentz boost. 

The final step is colour reconnection. In the old MI framework, Rick Field found a good 
agreement to CDF data if 90% of additional interactions produced two gluons with "nearest 
neighbour" colour connections [7j. In PYTHIA 8, with its more general MI framework, colour 
reconnection is performed by giving each system a probability to reconnect with a harder 
system 

2 

P = 7Z 2 ±R T^y P^c = RR * P±o, (8) 

{P±Rec + P±) 

where RR, ReconnectRange, is a user-tunable parameter and p± is the same parameter as 
in eq. ([3]) . The idea of colour reconnection can be motivated by noting that MI leads to many 
colour strings that will overlap in physical space. Moving from the limit of Nc — > oo to 
Nc = 3, it is perhaps not unreasonable to consider these strings to be connected differently 
due to a coincidence of colour, so as to reduce the total string length and thereby the 
potential energy. With the above probability for reconnection, it is easier to reconnect low 
P± systems, which can be viewed as them having a larger spatial extent such that they 
are more likely to overlap with other colour strings. Currently, however, given the lack of 
a firm theoretical basis, the need for colour reconnection has only been established within 
the context of specific models. 



3 Rescattering 

A process with a rescattering occurs when an outgoing state from one scattering is allowed 
to become the incoming state in another scattering. This is illustrated schematically in 
Figure [H where (a) shows two independent 2 — > 2 processes while (b) shows a rescattering 
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(a) 



(b) 



Figure 1: (a) Two 2 — > 2 scatterings, (b) a 2 — > 2 scattering followed by a rescattering 

process. An estimate for the size of such rescattering effects is given by Paver and Treleani 
[S] , where a factorised form is used for the double parton distribution, giving the probability 
of finding two partons of given x values inside an incoming hadron. Their results show that, 
at Tevatron energies, rescattering is expected to be a small effect when compared against 
the more dominant case of multiple disconnected scatterings. 

If we accept MI as real, however, then we should also allow rescatterings to take place. 
They would show up in the collective effects of MI, manifesting themselves as changes to 
multiplicity, p± and other distributions. After a retuning of p± and other model parame- 
ters, it is likely that their impact is significantly reduced, so we should therefore ask whether 
there are more direct ways in which rescattering may show up. Is there perhaps a region 
of low p± jets, where an event is not dominated by ISR/FSR, where this extra source of 
three-jet topologies will be visible? A further consideration is that such rescatterings will 
generate more p± in the perturbative region, which may overall mean it is possible to re- 
duce the amount of primordial k± and colour reconnections necessary to match data, as 
discussed in Section I2.4L 

3.1 Rescattering in Pythia 8 

If we begin with the typical case of small-angle t-channel gluon scattering, we can imagine 
that a combination of a scattered parton and a hadron remnant will closely match one of 
the incoming hadrons. In such a picture, we can write the complete PDF for a hadron as 



where the subscript u/5 is the unscattered/scattered component. That is, each time a 
scattering occurs, one parton is fixed to a specific x n value, while the remainder is still 
a continuous probability distribution. In such a picture, the momentum sum should still 
approximately obey 



Of course, in general, it is not possible to uniquely identify a scattered parton with one 
hadron, so an approximate prescription must be used instead, such as rapidity based. If we 



f(x,Q 2 ) -> frescaled(x,Q 2 ) + ^ 5 (x - X n ) = f u (x,Q 2 ) + fs(x,Q 2 ) 



(9) 



n 




(10) 
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Tevatron 


LHC 


Min Bias 


Jets 


Min Bias 


KclLsL) Jets 


Scatterings 


2.81 


5.11 


5.21 


12.20 


Single rescatterings 


0.37 


1.20 


0.93 


3.64 


Double rescatterings 


0.01 


0.03 


0.02 


0.11 



Table 1: Average number of scatterings, single rescatterings and double rescatterings in 
minimum bias and QCD jet events at Tevatron (y/s = 1.96 GeV, QCD jet p± m in = 20GeV) 
and LHC {y/s = 14.0 TeV, QCD jet p ±min = 50 GeV) energies 

consider the original MI probability given in eqs. ([T]) and , we can now generalise this to 
include the effects of rescattering 

dV M i ^ dV uu | dV nS | dV Su | dV S s ^ 
dp± dp± dp± dp± dp± ' 

where the uu component now represents the original MI probability, the u<5 and 5u compo- 
nents a single rescattering and the 55 component a double rescattering, where both incoming 
states to an interaction are previously scattered partons. 

Some indicative numbers are given in Table [U which shows the average number of 
scatterings and rescatterings for different types of event at Tevatron and LHC energies. The 
average distribution of such scatterings per event is also shown in Figure [2] for Tevatron 
minimum bias events. In the upper plot of dN/d(\ogp 2 ± ), the suppression of the cross 
section at small p\ is caused mainly by the regularisation outlined in eq. but is also 
affected by the scaling violation in the PDFs. Below p\ ~ 1 GeV 2 , the PDFs are frozen, 
giving rise to an abrupt change in slope. Normal scatterings dominate, but there is a 
clear contribution from single rescatterings. In the upper plot, it is not possible to see the 
effects of double rescattering, but this is (barely) visible in the ratio plot below. Given the 
overall small contribution from double rescatterings, we neglect these in the following. As 
previously predicted, rescattering is a small effect at larger p± scales, but, when evolving 
downwards, its relative importance grows as more and more partons are scattered out of 
the incoming hadrons and become available to rescatter. Note that here, we classify the 
original scattering and the rescattering by p^, but make no claims on the time ordering of 
the two. 

3.2 Mean p± vs Charged Multiplicity 

While a preliminary framework is in place which allows for hadronic final states, there are 
non-trivial recoil kinematics when considering the combination of rescattering, FSR and 
primordial k±. With the dipole-style recoil used in the parton showers, a final-state ra- 
diating parton will usually shuffle momenta with its nearest colour neighbour. Without 
rescattering, colour dipoles are not spanned between systems, and individual systems will 
locally conserve momentum. With rescattering enabled, you instead have the possibility 
of colour dipoles spanning different scattering systems and therefore the possibility of an 
individual system no longer locally conserving momentum. When primordial k± is now 
added through a Lorentz boost, these local momentum imbalances can lead to global mo- 
mentum non-conservation. In order to proceed and be able to take an initial look at the 
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Figure 2: Average distribution of scatterings, single rescatterings and double rescatterings 
per event (y/s = 1.96 GeV, minimum bias). Double rescattering is not visible at this scale 
in the diV/d(logp5_) plot, but is visible in the ratio 

effects of rescattering on colour reconnection, a temporary solution of deferring FSR until 
after primordial k± is added has been used, as is done in Pythia 6.4. 

We begin by studying the mean p± vs charged multiplicity distribution, (p±)(n c h), from 
Pythia 6.418 (Tune A) and Pythia 8.114 (default settings), compared to the CDF Run 
II data (\rj\ < 1 and p± > 0.4 GeV/c) [§]. For each run, the p± parameter of the MI 
framework is tuned so that the mean number of charged particles in the central region 
is maintained at the Tune A value. This is shown in Figure El where we can see that 
Pythia 6, using virtuality-ordered showers and the old MI framework, does a reasonable 
job of describing the data. Pythia 8 does not currently have a full tune to data, but does 
qualitatively reproduce the shape of the data when colour reconnection is turned on, up to 
an overall normalisation shift. It is clear that without colour reconnection, the slope of the 
curve is much too shallow and unlikely to describe the data, even given an overall shift. 
The same results with deferred FSR are also shown; the slope is marginally steeper, but 
still in the same region as without deferred FSR. 

FigureHJnow shows the results when rescattering is enabled. Starting without any colour 
reconnection, we see that when rescattering is turned on, there is a rise in the mean p±, 
but also that this is in no way a large gain. This is also the case when colour reconnection 
is turned on and tuned such that the curve qualitatively matches the shape of the Run II 
data. The amount of colour reconnection used is given in the form RR*p± , as described in 
eq. (jSj). That a rise in the mean p±_ is there with rescattering, but small, is something that 
was observed already in an early toy model study. Now, when the full generation framework 
is almost there, it is clear that rescattering is not the answer to the colour reconnection 
problem. Other potential effects of rescattering remain to be studied. 
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Figure 3: Mean p± vs Charged Multiplicity, \rj\ < 1 and p±_ > 0.4 GeV/c, CDF Run II 
data against Pythia 6.418 (Tune A) and Pythia 8.114 (default settings) with and without 
deferred FSR 
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Figure 4: Mean p± vs Charged Multiplicity, Pythia 8.114 (deferred FSR), effects of 
rescattering 
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4 Enhanced Screening 



The idea of enhanced screening came from the modelling of initial states using dipoles in 
transverse space [TU]. A model using an extended Mueller dipole formalism has recently 
been used to describe the total and diffractive cross sections in pp and 7*p collisions and the 
elastic cross section in pp scattering [IT]. In such a picture, initial-state dipoles are evolved 
forwards in rapidity, before two such incoming states are collided. In the model, as the 
evolution proceeds, the number of dipoles with small transverse extent grows faster than 
that of large dipoles. The dipole size, r, determines the screening length, which appears 
in the interaction cross section as a p± cutoff, p± ~ 1/r. Smaller dipoles imply a larger 
effective cutoff, and an enhanced amount of screening. A rough calculation shows that this 
screening effect is expected to grow as the square root of the number of dipoles. 

To model this in Pythia, we consider the p± parameter of the MI framework that 
encapsulates colour screening, as given in eq. ([3]). By scaling this value by an amount that 
grows as the amount of initial-state activity grows, this enhanced screening effect can be 
mimicked. Such a change can be achieved by adjusting the weighting of the cross section 
according to 

^ a 2 s {p 2 ±0 + pl) _^ a 2 s {P 2 ±o + Pl) (12) 

dpi (pio + pi) 2 (npio+pi) 2 ' 

where n takes a different meaning for two different scenarios. With the first scenario, ESI, 
n is set equal to the number of multiple interactions that have taken place in an event 
(including the current one). In the second, ES2, n is set equal to the number of MI+ISR 
interactions that have taken place in an event. 

4.1 Mean p± vs Charged Multiplicity 

We again study the (p±)(n c h) distribution, this time with the enhanced screening ansatz. 
The results are given in Figure El Looking at the curves without colour reconnection, it is 
immediately apparent that both scenarios give a dramatic rise in the mean p±, although 
not quite enough to explain data on their own. With colour reconnection now enabled 
and tuned, again so that the curves qualitatively match the shape of the Run II data, it is 
possible to noticeably reduce the amount of reconnection needed. With colour reconnection 
at these levels, there is still perhaps an uncomfortably large number of systems being 
reconnected, but the results are definitely encouraging. There are many more areas to 
study in relation to enhanced screening, but from these initial results, it is worth checking 
if it may play a role in reducing colour reconnections to a more comfortable level. 

5 Conclusions 

Pythia 8, the C++ rewrite of the Pythia event generator has now been released. It 
has been written with a focus on Tevatron and LHC applications, something that is ev- 
ident given the sophisticated MI model present in the program. The original MI model, 
introduced in the early versions of Pythia, has been well proven when compared to ex- 
perimental data. The new Pythia 8 MI framework, based on this original model, now 
generalises the physics processes available, as well as adding entirely new features. 
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Figure 5: Mean p± vs Charged Multiplicity, Pythia 8.114, effects of the enhanced screen- 
ing ansatz 

We have also taken an early look at rescattering and enhanced screening, two new 
ideas for modifying the physics inside the MI framework. There is currently a preliminary 
framework for rescattering, although fully interleaved ISR, FSR and MI is still to come. It 
appears, at this early stage, that rescattering is not the answer to the colour reconnection 
problem, but there is still much more to investigate, such as three-jet multiplicities and other 
collective effects. The idea of enhanced screening leads to a simple ansatz that gives large 
changes when looking at the (p±)(n c f l ) distribution. Again, there are still many questions 
to be asked, including how this modification affects other distributions. 
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